syms a b c d
D = [a b c ; a b d; d c b; c b a];
d = size(D)
r = size(D, 1)
